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ABSTRACT 

We derive explicit Lorentz-invariant solution of the Einstein and null geodesic equations for 
data processing of the time delay and ranging experiments in gravitational field of moving 
gravitating bodies of the solar system - the Sun and major planets. We discuss general- 
relativistic interpretation of these experiments and the limitations imposed by motion of the 
massive bodies on measurement of the parameters Yppn, Pppn and 8ppN of the parameterized 
post-Newtonian formalism. 
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1 INTRODUCTION 



Theoretical speculations beyond the Standard Model suggest that gravity must be naturally accompanied by a partner - one or more scalar 

fields , which contribute to the hybrid metric of space-time through a system of equations of a scalar-tensor gravity theory jPamour & Esposito-Farese 
19921) . Such scalar partners generically arise in all extra-dimensional t heories, and notab ly in string theory. Scalar fields play also an impor- 
tant role in modern cosmological scenarios with the inflationary stage ( iMukhanovluOOSr) . Therefore, unambiguous experimental verification 
of existence of the scalar fields is among primary goals of gravitational physics. 

Phenomenological presence of the scalar field in the metric tensor is parameterized by three parameters - Yppn, Pppn and Sppfj - of the 
parameterized post-Newton ian (PPN) formalism. These parameters enter the metric tensor of a s tatic and spherically-symmetric gravitating 
body in the following form JPamour & Esposito-Faresell 19961 : limmbergll 19921 : IwilJ 19931 . [200 lb 
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where we have used the isotropic coordinates X 
PPN parameters fppN = ypp^ — 1, Pppn = Pppn 



= (cT,X), R = \X\, and denoted deviation from general relativity with the c omparative 
1; SppN = dp pN— 1- Parameter 5ppN generalizes the standard PPN formalism I Willi 1993 ) 



to the second post-Newtonian approximation iBrumberg 1992 ). One notices that 5p/>/v is actually related to pppff and yppN in a generic 
scalar-tensor theory of gravity i Damour & Esposito-Faresel 19961) . In particular, this theory predicts that jSpp/v cancels in the combination 
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^PppN + -^/^SppN entering equation. (4) of tlie present paper, which should depend, theoretically, only on yppN and its square. Thus, high- 
precision missions will have a very clean access to Yppn- However, we prefer to keep combination —pppf^ + 3/4Spp!\i explicitly in our 
equations in order to separate parametrization of the second post-Newtonian effects associated with SppN from the linearized Shapiro time 
delay, which is parameterized by Yppn alone. Moreover, parameter SppN is independent from fippN and Yppn in vector-tensor theories of 
gravity i Deng. Xie & Huan^boogh . In general relativity, pppN = Yppn = SppN = 0. 

The best experimental boun d on Yppn = (2.1 ±2.3) x 10~ ^ has been obtained (under a certain implicit assumption I lKopeikin et al 



2007^ ) in the Cassini experiment l IBertotti. less & Tortorall 2003^ . Limits on the parameter Pppn depend on the precision in measuring fpPN, 
and are derived from a linear combination IfppN — &PP N < 3 x 10~^ by observing the Mer cury's perihelion shift, and from AfippN — Yppn = 
(4.5 ± 4.5) X 10^^ imposed by the lunar laser ranging ( Williams. Turvshev & Bogg3l2004 ). Parameter 5ppN has not yet been measured. 

The most precise measurement of Yppn and 8ppN can be achieved in near-future gravitational experiments with light propagating in the 
field of the Sun or a major planet. Post-post-Newtonian e quation of the relativistic time delay in a static gravitational field is obtained from 
the m etric (TJ, l[2). It was derived by a number of authors i Richter & Matzneilll982 , 1983 ; Brumber3ll992 ; Tevssandier & Le Poncin-Lafitte 
200 sh and reads (in the isotropic coordinates) as follows 
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where Ti and T2 are coordinate times of emission and observation of photon, R - 
emission, Xi, and observation, X2, of the photon, and 
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is the extra time delay caused by the gravitational field, Ni = X\/Ri and N2 = X2/R2 are the unit vectors directed outward of the 
gravitating body, R\ = \X\\, R2 = IX2I are radial distances to the points of emission and observation respectively. 

The Sun and planets are not at rest in the solar system because they are moving with respect to the barycenter of the solar system as 
well as with respect to observer. Motion of the light-ray deflecting body (the Sun, a major planet) affects propagation of light bringing the 
post-Newtonian corrections of the order of (GM/c^){v/c), [GM / c')[v / c)^ , etc. to equation where v is a characteristic speed of the 
massive body with respect to a reference frame used for data processing, which can be chosen as either the barycentric frame of the solar 
system or the geocentric frame of observer. These motion-induced post-Newtonian corrections to the static time delay Ar correlate with the 
PPN parameters making their observed numerical value biased. Therefore, it is important to disentangle the genuine effects associated with 
the presence of the scalar field from the special-relativistic effects in equation l|4ll imparted by the motion of the bodies. 

This problem has not been addressed until recently becau se the accuracy of astronomical observa tions was not high enough. How ever, 
VLBI measurement of the null-cone gravity-retardation effect 1 KopeikiJboOl : Fomalont & KoDeikinll2003 ; ICopeiMn 20(34; F omalont et al.l 
2OO9I) and frequency- shift measurement of Yppn in the Cassini experiment jAnderson. Lau & Giampieril2004l : lBertotti. Iess^& Tortoral 2003h 
made it evident that modern technology has achieved the level at which relativistic effects c aused by the dependence of t he gravitational 
field on time can be no lo nger ignored. Futur e gravitational light-ray deflection experim ents 1 Kopeikin & Mashhoori 2002), radio ranging 
BepiC olombo experiment l IMilani et al.ll2002h . laser ranging experiments ASTROD ( INiluOOTl) and LATOR l ITurvshev, Shao & Nordtvedt 
20040 will definitely reach the precision in measuring Yppn-, Pppn and SppN that is comparable with the post- Newtonian corrections to th e 
static time delay and to the deflection angle caused by the motion of the massive bodies in the solar system 1 Plowman & Hellings 20061) . 
Therefore, it is worthwhile to undertake a scrutiny theoretical study of the time-dependent relativistic corrections to the static Shapiro time 
delay. 

In this paper we focus on deriving two apparently different forms of the Lorentz invariant solution of the light ray equations (see 
equations (31) and (43)) in the linearized (with respect to the universal gravitational constant G) approximation of general relativity by making 
use of the technique of the Lienard-Wiechert potentials 1 Kopeikin & Schafeilll999l) and algebraic transformations of th e retarded quantities. 
In particular, equation (43) of the present paper significantly generalizes the result of Bertotti. Ashbv & lessl 1 20081) for the gravitational 
time delay. We expand this retarded-time solution in the post-Newtonian series in three various ways (see equations (63), (68) and (86) 
below) and analyze the impact of the velocity-dependent corrections on measuring values of the PPN parameters in the gravitational time- 
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delay experiments. Section VIII discusses a correspondence between the Lorentz symmetry group for gravity and light as revealed by the 
time delay experiments. Section IX gives a justification that the ODP code of NASA must be revamped for doing adequate processing of 
high-precise data in ranging gravitational experiments. 



2 NOTATIONS 

In what follows the Greek indices a,j3, ... run from to 3, the Roman indices /, /, ... run from 1 to 3, repeated Greek indices mean Einstein's 
summation from to 3, and bold letters a = (a} ,0^ ,a}),b = {b^ ,b^,b^), etc. denote spatial (3-dimensional) vectors. A dot between two 
spatial vectors, for example a • 6 = fl'i' + a^h^ + a^b^, means the Euclidean dot product, and the cross between two vectors, for example 
axb, means the Euclidean cross product. We also use a shorthand notation for partial derivatives da = d/dx". Greek indices are raised and 
lowered with full metric gg,p- The Minkowski (flat) space-time metric rjap = diag(— 1, + 1, + 1, + 1). This metric is used to rise and lower 
indices of the unperturbed wave vector of light, and the gravitational perturbation h^p ■ 



3 THE LIENARD-WIECHERT GRAVITATIONAL POTENTIALS 

We introduce the post-Minkowskian decomposition of the metric tensor 

Sap = ^ap +^>ap > 



(5) 



where h^p is the post-Minkowskian perturbation of the Minkowski metric tensor Tj„p . We impose the harmonic gauge condition jMisner. Thome & Wheelei 
1973h on the metric tensor 



dah"P-^dPh\=0. 

In arbitrary harmonic coordinates x" = [ct, x), and in the first post-Minkowskian approximation the Einstein equations read 



c2 dt'- 
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where T^^ is the stress-energy tensor of a light-ray deflecting body. In linearized approximation this tensor is given by the following equation 
(f , x) = Muf^u^^Jl - i325(^> (a; - z{t)) , (8) 

where M is the (constant) rest mass of the body, z (r) is time-dependent spatial coordinate of the body, l3 = c^^ dz/dt is velocity of the body 
normalized to the fundamental speed c. 



-1/2 



-1/2 



(9) 



is the four-velocity of the body normalized such that iiaii'^ = — 1, and S'-'' (x) is the 3-dimensional Dirac's delta-function. We have neglected 
^— g in equation ^ because in the linearized approximation \/—g = I + 0{G), and the quadratic terms proportional to are irrelevant 
in T^^ since they will give time-dependent terms of the second post-Minkowskian orde r of magnitude, which are curre ntly negligible for 
measurement in the solar system. For the same reason, we do not use the metric derived by Blanchet. Fave & Ponsoll i I998h as it goes beyond 
the approximation used in the present paper. We have also used a standard notation f3 for the dimensionless velocity of the body. This notation 
should not be confused with the PPN parameter fippN- 

Because the Einstein equations ^ are linear, we can consider their solution as a linear superposition of the solutions for each body. It 
allows us to focus on the relativistic effects ca used by one bod y (the Sun, planet) only. Solving Einstein's equations Q b y making use of 
the retarded Lienard-Wiec hert tensor potentials Bel et aljl981 ), one obtains the post-Minkowski metric tensor perturbation teel etal.ll98ll 
Kopeikin & Schafeilll999h 
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where 
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In equation ilOi all time-dependent quantities are taken at retarded time j defined by the null cone equation ( I13t given below, u" 



dz"{s)/ds is its four- velocity, with being a retarded time (see below), l3{s) 



dz{s)/ds is body's coordinate velocity normalized 



to the fundamental speed c. Notice that the metric tensor perturbation JlOb is valid for accelerated m otion of the gravitating body as well, 
and is not restricted by the approximation of a body moving on a straight line (see lBel et alj jl98lh for more detail). In other words, the 



four-velocity in equation jlOt is not a constant, taken at one, particular event on the world line of the body. 

Because we solved the Einstein equations Q in terms of the retarded Lienard-Wiechert potential s, the distance p 
body's worldline z"(.s) = {cs,z{s)), and the four- velocity are all functions of the retarded time s i Bel et aljl98lh . The retarded time J 
is found in the first post-Minkowski approximation as a solution of the null cone equation 

x^-z^.))(x^-z^(.))=0, 



the 
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where the constant c in equation l ll4t denotes the fundamental speed in the Minkowski space-time, which physical meaning in equation ( 114b 
is the speed of propagation of gravity as it originates from the gravity field equations ([7}. It is important to notice that equation l ll4t is a 
complicated function of t he retarded time s = s{t,x), which has an analytic solution only in case of a uniform motion of the gravitating 
body along a straight line jKoDeikiij|2004l) . Geometrically, equation l lI4t connects the point of observation x and the retarded position of the 
gravitating body z(s) by a null characteristic of the linearized Einstein field equations (O. Radio waves (light) are also propagating along a 
null characteristic connecting the observer and the radio emitter. However, the null characteristic of the linearized Einstein equations il4\ is 
well separated on the space-time manifold (and in the sky) from the null characteristic associated with the propagation of the radio wave in 
any kind of ranging and time-delay experiments. Hence, they should not be confus ed in relativistic experiments involving l ight propagation 
in the field of a moving gravitating body, which gravitational field depends on time ( lwili200ll : lKopeikin & Fomalonllbooeh . 

All components of the time-dependent gravitational field (the metric tensor perturbation h^p) of the solar system bodies interact with 
radio (light) waves moving from a radio (light) source to the Earth, and perturb each element of the phase of electromagnetic wave with the 
retardation given by equation J14b . The use of the retarded Lienard-Wie chert gravitational potentials, rather than the advanced potentials, 
is consistent with the principle of causality i Kopeikin & Fomalonjl 2007 ), and the observation of the orbital decay of the re lativistic binary 
pulsar B1913+16 caused by the emission of gravitational radiation, according to general relativity i Weisberg & Tavlo3l2005 ). 



4 THE ELECTROMAGNETIC PHASE 



Any ranging or time delay experiment measures the phase y of an electromagnetic wave coming from a spacecraft or a radio (light) source 
outside of the solar system. The phase is a scalar function bei ng invariant with respect to coord i nate transformations. It i s determined in the 
approximation of geometric optics from the eikonal equation I Misner, Thome & Wheeleilll973l ; Landau & Lifshitj|l97I ) 



8l'''df,xi/dvV = 0, (15) 

where g^'^ = T)^'^ — h^^ . T he eikonal equation dlSt is a direct consequence of Maxwell's equations i Misner. Thome & Wheelei I973I : 
Kopeikin & Mashhooijliooj) and its solution describes localization of the front of an electromagnetic wave propagating on a curved space- 
time manifold, which geometric properties are defined by the metric tensor ([S}, l llOt that is a solution of the Einstein equations. We emphasize 
that the electromagnetic wave in equation l ll5t has no back-action on the properties of the metric tensor g^v, and does not change the curvature 
of the space-time caused by the presence of the gravitating body. Thus, experimental studying of the propagation of electromagnetic wave 
allows us to measure the important properties of the background gravitational field and space-time manifold. 
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Let us introduce a co- vector of the electromagnetic wave, Ka = da iff- Let A be an affine parameter along a light ray being orthogonal 
to the electromagnetic wave front ij/. Vector K" = dx" /dX = g"^dp is tange nt to the light ray. Equati on l ll5t expresses a simple fact that 
vector A"" is null, that is gjivK^K'^ — 0. Thus, the light rays are null geodesies i Landau & Lif shit3l 1 97 1 ) defined by equation 

'-^^{dag^.K^K^. (16) 

The eikonal equation J15t and light-ray equation l ll6b have equivalent physical content in general relativity since equation ( I15t is a first 
integral of equation J16t . 

Regarding propagation of electromagnetic wave, it is more straightforward to find solution of equation dlSb . To this end, we expand the 
eikonal i// in the post-Minkowskian series with respect to the universal gravitational constant G assuming that the unperturbed solution of 
equation l ll5t is a plane electromagnetic wave (that is, the parallax of the radio source is neglected). The expansion reads 

V/=i//o + ^M" + <P(x")1 + 0(g2), (17) 

where xffQ is a constant of integration, k" = (l,fc) is a constant null vector directed along the trajectory of propagation of the unperturbed 
electromagnetic wave such that rj^vk^k'^ = 0, V is the constant frequency of the unperturbed electromagnetic wave, and (p is the first post- 
Minkowskian perturbation of the eikonal, which is Lorentz-invariant. Substituting expansions (jSj, OZJ to equation \l5h and leaving only 
terms of order G, one obtains an ordinary differential equation for the post-Minkowskian perturbation of the eikonal. 



V/^,„,, = £i^l!l^, (18) 



l[\n{-kap")] , (19) 



- = -h nakp = — — 

which can be also obtained as a first integral of the null geodesic equation ( I16t . Equation il8b can be readily integrated if one employs an 
exact relationship 

dX ds 1 

Pr kap" kaifi 

which makes the integration straightforward. Indeed, if the body's acceleration is neglected, a plane-wave solution of equation dlSb is 

^ 2GMV (^^^,a)i„(_^^pa) ^ (20) 

where all quantities in the right side are taken at the retarded instant of time s in compliance with the null cone equation il4i . One notices that 
the time of the closest approach of the light ray to the moving body does not play any role in calculation o f the gravitational perturbation 
of the electromagnetic phase. The time is a good approximation of the retarded time s jKopeikin 2001 ), and ca n be used in practical 
calculations of light propagation in the gravitational field of moving bodies jKlioner & Kopeikiij|l992l : lKlionerll2003bh . However, it does not 
properly reflect the Lorentz-invariant nature of the gravitational time delay and makes its post-Newtonian expansion looking more entangled 
and complicated. Further discussion of this issue is given in section IT!2l 

One can easily check that equation l l20b is a particular solution of equation l ll5b . Indeed, observing that 

daP" =5^-u^'daS, (21) 

one obtains from the null cone equation l ll3b 

das = -^. (22) 
PR 

Differentiation of equation ( |20b using equations ( 121b and ( 122b shows that equation dlSb is satisfied. 

Equation l l20b for the electromagnetic phase is clearly Lorentz-invariant and valid in an arbitrary coordinate system. It tells us that a 
massive body (the Sun, planet) interacts with the electromagnetic wave by means of its gravitational field, which originates at the retarded 
position z{s) of the body and propagates on the hypersurface of null cone ( 114b . The gravitational field perturbs the phase front of the 
electromagnetic wave at the field point x" regardless of the direction of motion of the incoming photon or the magnitude of its impact 
parameter with respect to the body. This consideration indicates a remarkable experimental opportunity to observe the retardation effect of 
the gr avitational field by measuri ng the shape of the ranging (Shapiro) time delay and comparing it with the JPL ephemeris position of the 
body Jstandish & Williamsl2006l) obtained independently from direct radio/optical observations of the body, conducted in preceding epochs. 
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This idea was executed in VLBI experiment with Jupiter jKopeikijlioOll uFomalont & Kopeikinllioosi) . Next section explains the null-cone 
relationship between the characteristics of the Maxwell and Einstein equations. 



5 THE RANGING TIME DELAY 

The Lorentz-invariant, general-relativistic time delay equation, generalizing the static Shapiro delay l lShapiroll964l) . can be obtained directly 
from equation J20b . We consider a ranging time-delay experiment in which an electromagnetic wave (a photon) is emitted at the event with 
4-dimensional coordinates x" = (c/j , a; i ) , passes near the moving gravitating body, and is received by observer at the event with coordinates 
x" = {ct2,X2)- In the most general case, the emitter and observer can move, which means that coordinates xi and X2 must be understood as 
functions depending on time ti and t2 respectively, that is xi = x{ti) and X2 = x{t2), where x{t) is a spatial coordinate of the photon taken 
at time t. The gravitating body is also moving during the time of propagation of the electromagnetic wave from the emitter to the observer. 
In the approximation of a uniform and rectilinear motion, which is sufficient for our purpose, spatial coordinate of the body is given by a 
straight line 



z{t) = ZQ+Vt . 



(23) 



where zq is position of the body taken time f = 0. One notices that the spatial coordinate of the body entering the Lienard-Wiechert solution 
of the gravity field equations depends on the retarded time s. It means that the time argument t in equation ( I23t must be replaced with the 
retarded time s without changing the form of this equation. In other words. 



: ZO + VS , 



(24) 



where the retarded time s is given by the solution of the gravity null cone equation il4t . In case of a rectilinear and uniform motion of the 
gravitating body 



R 13 



(25) 



c(l-/32) 

and R = X ~ z{t) with z{t) defined in equation ( I23t . 

The unperturbed spatial components (^') = fc of the wave vector fc" are expressed in terms of the coordinates of the emitting and 
observing points 

X2-XX 



k 



(26) 



\X2~Xl\ 

This vector is constant for a single passage of the electromagnetic wave from the emitter to the observer. However, in case when the emitter 
and/or observer are in motion, the direction of vector k will change as time prog resses. This remark is important for calculation of the 
Doppler shift of frequency, where one has to take the time derivative of the vector k i Kopeikin & Schafe3l 19991 : Kopeikin et alj[^007 ) 



The perturbed wave vector, K" = dx" /dX, is obtained from the eikonal equation ilOi by making use of identification K" = dxfr/dx", 
which is a consequence of the Hamiltonian theory of light rays and can be used for further integration in order to determine the tra- 
jectory of propagation of th e electromagnetic wave in the curved space-time. The explicit integration has been performed in paper by 
Kopeikin & Fomalonll 1 20061) and could be used for calculation of the ranging time delay. However, in the present paper we shall rely upon a 
different method. 

We note that the phase xjf of the el ectromagnetic wave, emitted at the point x" = (cf|,a:i) and received at the point x^ = (cf7,a:?), 
remains constant along the wave's path i Misner. Thome & Wheelei 1973 : Landau & Lifshitz 1971 : Kopeikin & MashhoonI 2002 ). Indeed, 
since X is an affine parameter along the path, one has for the phase's derivative 

dxif dw dx" „ 

which means that y/(x"(A)) =const., in accordance with our assertion. Equating two values of the phase at the point of emission of the 
electromagnetic wave, x", and at the point of its receptions, x", and separating time from space coordinates, one obtains from equations (TT), 
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where the retarded distances p" 

1, , 



-z"(.^2),pr 



= — z" (ii ), and the retarded times S2, s\ are defined by the null cone equations 
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which are inferred from equation l ll4t . Expanding all Lorentz-invariant scalar products, and replacing relationship J26t in equation ( 
the ranging delay 



At 



1, 

-F2 

c 



xi\+At , 
2GM l-k /3 



In 



P2-kp2 



Pi-kpi 



(28) 

(29) 

(30) 
I yields 

(31) 
(32) 



where the retarded, null-cone distances p2 = X2 — z{s2), p\=x\ — z(s\), P2 = |p?.|. Pi = |pi I- 

Lorentz-invariant expression for ranging delay l l32t was derived first by iKopeikin & Schafer J 19991) by s olving equations for light 
geodesies in the gravitational field of moving bodies with the Lienard-Wiechert gravitational potentials. Later on, iKlionej J20033) obtained 
this expression by making use of the Lorentz transformation of the Shapiro time delay (which is equivalent to a simultaneous transformation 
of the solutions of both the Einstein and Maxwell equations) from a static frame of the body to a moving frame of observer. Notice that 
in general relativity equation i3l\ describes a hypersurface of the null cone along which both electromagnetic and gravitational field are 
propagating. Electromagnetic characteristic of the null cone is given by the null vector k of the photon, while the null characteristic of the 
gravity field enters the time delay equation l l32t in the form of the retarded time .?, which is the time argument of the coordinate z of the 
moving body under consideration. 

In the present paper we derive another useful form of the Lorentz-invariant expression for the ranging delay, which can be directly 
compared with and generalizes the approximate ranging delay formula currently used in the NASA Orbit Determination Program (ODP). 
This derivation comes about from the following exact relationships 



P2 - fe • P2 = 

Pi-fcpi = 
where r = |r|, r : 
r" = rk" = (r,r) 



\pi-zis2) + z{si)\^-{r-p2f 



2r 

\p2 + z{s2)-z{si) 



-Pi) 



2r 



- X2 — x\,s,Q that 



(33) 
(34) 

(35) 



is a null vector in the flat space-time connecting coordinates of the point of emission and reception of the electromagnetic wave: T]^^ r^^ =0. 
Because the gravitating body moves uniformly with constant speed v, its coordinate z{s) is not constant and can be expanded as follows (see 
equation l l24t ) 



z{s2) = z{si)+v{s2-si) , 

where the time interval ^2 — -^l can be expressed in terms of the null-cone distances by making use of the retarded time equations i29\ , 
and the ranging equation i3H . One has, 

1 



.52 -.51 



(i2-?2) + (f2-'l) + ('l -.5l) = -(r + pi -p2) + 0(c ^) 



(36) 
(37) 



Plugging equation l l37t to l l36t , and replacing it in equations i33\ . i34\ allows us to transform the ranging time delay logarithm to the following 
form 

P2 + Pi + 2(P2 • /3) - (r + pi - P2) 



In 



'P2 


-fe 


P2 


= -In 


.Pi 


-fe 


Pi. 





P2 + Pi - 2(pi ■ /3) + /32 (r + Pi - P2) 



(38) 
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Let us now make use of definition Jilt of the Lorentz-invariant distances 

P2« = -"aP2 = ^JYZT^ ' ^^^^ 



PIR = -«„p« = ^li^^^^ . (40) 



Pi 

Tedious but straightforward calculations reveal that 



p2 + Pl+'--2(p2-/3)-|8n'- + Pl-P2) = \/l-i32(P2fi + Pl«-'-fca«") , (41) 



P2+Pi-r-2(pi-/3)+j8-(r + pi-p2) = \/l-/32 (pjx + PxR + rka ) . (42) 
These equations taken along with equation ( 135b allows us to reduce the time delay logarithm in equation l l38b to another Lorentz-invariant 

form 

/ P2-k-P2\ P2K+Pm-Pn\ (43^ 



VP2fi + PlR+Pl2 

where the ranging distance pi2 = rkaii" = uaf"' is invariant with respect to the Lorentz transformation. It represents contraction of the null 
vector defined in equation ( I35t with four-velocity h" of the gravitating body. The null vector determines (unperturbed) propagation 
of the electromagnetic signal. Distances pijj, P2R are defined in equations l |39t . l |40t . and they also represent contraction of the null vectors 
p", p" with four-velocity of the gravitating body. However, contrary to vector r", vectors pf, p" describe the null characteristics of the 
gravitational field. 

Accounting for equation ( 143 1 the Lorentz-invariant expression for the time delay assumes the following form 
^^^2GMl^Jp,, + p,,-p,,\ ^^^^ 



Vl-jS^ VP2R + PiR+P12 

This equation is apparen tly Lorentz-invariant, valid for any value of the velocity of the light-ray deflecting body, and essentially generalizes 
the result of the paper by Bertotti. Ashbv & Ies3 1 200 sh . 



6 POST-NEWTONIAN EXPANSION OF THE RANGING DELAY 

Let us introduce an auxiliary vectors dBel et al .Il98lh 



^ d"P2R ■ 



P2_ 
P2R 



d"Pl 



PT 
PlR 



(45) 



Vectors p" and p" are null as defined by the (gravity-field) null cone equations ( 129b . ( 130b . The four- velocity of the body, w", is a time-like 
vector. Halt" = — 1. The difference between the null and time-like vector yields the space-like vectors n", ii", because nian" = nia'T-^ = +1- 
The post-Newtonian expansion of z"(i2) around time f2, and the post-Newtonian expansion of z"{si) around time ti are obtained by 
making use of a Taylor expansion. Omitting acceleration, one gets 



P2" 

Pf 
and 

P2 
Pi 



'^2 ' 



{s2-t2) 



ds 



P2R 
PlR 



as 



- Up 



■r^+P2i'" 



r'^+piu'' 



p 



(46) 
(47) 

(48) 
(49) 
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where the retarded time equations i29\ , i30l have been used to replace time intervals ^2 — ^2 ^nd ^ i — f i . We have also introduced in previous 
equations the pure spatial vectors 

= 4^z"{t2) = {4 = 0,r'2=xi,-z'{t2)} , (50) 
rf = ^f-z"(t,) = {r'l = 0,r\=x\-zHh)} , (51) 

which are lying on the hypersurface of constant time to and ti respectively. 
Substituting equations (I46t- ([49} into equation J45t reveals that 

«2P2R = '-2 '•2). (52) 

nfpiR = r« + H«(«/jrf). (53) 
Taking into account that n" and are space-like unit vectors, one has 



P2R = \l'-axr<^ + {^ar1Y = \l'^-^^^, (54) 



PlR = ^Jraorf + {uarf) = ^ ''' ^'^''^P^' . (55) 
We further notice that, if acceleration is neglected, 

P12 = ,^ '^^^ rn , (56) 

vi-p 

where the unit vector 
k-f3 



\k-f3\ ' 
the relative distance 



(57) 



ri2 = \r2-ri\, (58) 
and 

r2 = X2-z{t2), (59) 

ri = xi-z{ti), (60) 

are spatial distances from the observer to the body and from the emitter to the body taken respectively at the time of reception and that of 
emission of the electromagnetic wave. It is worth observing that the post-Newtonian expansion of the Euclidean dot product k ■ a does not 
have a term, which is linear with respect to velocity 

fe-fT = l-i(fcx/3)2 + 0(jS3). (61) 
This expansion yields 

Pi2 = n2 + 0(/3-), (62) 

that is the distance r\2 is a Lorentz-invariant function up to the second post-Newtonian corrections of the order of jS". This justifies the 
replacement of the heliocentric c oordinates of the massive bodies of the solar system to their barycentric counterparts introduced by Moyer 
in the ODP manual 1 MoveillioOS ) ad hoc (see section|9]of the present paper for further details). 

After preceding preparations, we are ready to write down the post-Newtonian expansion for the ranging time delay. We would like 
to emphasize that the post-Newtonian expansion of the ranging delay is not unique and can be represented in several different forms, 
which are physically and computationally equivalent. However, this non-uniqueness complicates things and has been debated in papers 
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I Bertotti. Ashbv & lesslbooi : Kopeikinl 20091) regarding the nature of the relativistic time delay effects associated with motion of the grav- 
itating body. In what follows, we derive all possible forms of the post-Newtonian expansion of the ranging delay demonstrating that the 
relativistic effects associated with the motion of the ligh t-ray deflecting body are induced by the grav itomagnetic field arising due to the 
translational motion of the body with respect to observer i Kopeikin & Mashhoonll2002 ; lKopeikiij|2004h . 
First of all, substituting equations ( I39l l. ( I40t to ( 143 1 casts the ranging delay B2b in the following form 



2GM 1 - fc • /3 



In 



rl - (/3 X roY + Jr\-{(ix n )2 + (fc • CT)rn 



^rl - (/3 X r2)2 + - (/3 x rj )2 - (fc ■ a)rn 



(63) 



which is the most convenient for making its explicit post-Newtonian expansion with respect to the ratio of j6 = v/c. Neglecting terms of the 
order of one has 



+ 



l-fc-/3 + i/32 

GM rn (ni x /3) 

3 



2GM, 



■In 



(64) 



c- rir2 
where the unit vectors n i 



2n + (n2X/3)2r2-(fcx/3)2(ri+r2) ^ g ( GM 

1+ni ■n2 V ^'^ J 



ri /ri, no = t'2/''2 with ri, r2 being defined in equations l |59b , J60t (see Fig.lTJ. 
Velocity-dependent corrections appear in this expression explicitly as the terms depending on P = v / c, and implicitly in the argument 
of the logarithm, which depends on two positions of the body taken at times t\ and t2, that is z{t2) = z{t\) + v(t2 — ti) = z{ti) + /3r so that 
ro and ri2 are not independent of rj. We discuss the impact of the velocity-dependent terms on measured values of the PPN parameters in 
the next section. 

It is also instructive to derive the time delay equation in the linearized form as it is given in 



Will 



20011) . We make use of equations 



( I45l l- J55b to get the post-Newtonian expansion of functions entering the argument of the logarithm in the ranging delay ( 
kap" 



kaPi 



kar^ + {kau") 



karf + ikaii" 



r2pr2 + ("P'2 



"P'l + 



IP'^I 



'P'i 



(65) 



(66) 



Explicit expansion of these equations with respect to the powers of the velocity-tracking parameter jS = v/c brings about the following result 
r2 - fc • r2 +/3 ■ r2 - r2 (fc • /3) + O f j8 



P2-k p2 = r2-fc-r2+/3-r2-r2(fc-/3) + 0(j8^) , (67) 

Pi-kpi = n-k-r^+f3-ry-n{k-fi) + 0[fi^y (68) 

Applying these expansions to the argu ment of log arithm in the ranging delay | |32I ( yields the first term in the post-Newtonian expansion of 
the ranging delay in the form given in i Wili20oi ) 



Af : 



(1- 



,2GM 
-fc-/3)— ^In 



^2 - (T • r2 



ri -cr-ri 



O 



2GMp- 



(69) 



where the unit vector 

(T = fc-fex(/3xfe) + 0(iS-), (70) 

is the same as that defined by equation iSli . 

The explicit post-Newtonian dependence of the time delay on velocity of the gravitating body v enters the argument of the logarithm 
in the form of equation ( 170b . which looks like the aberration of light for the unit vector fc. However, equation ( 169b approximates the exact 
time delay equation ( 131b . which demonstrates that the argument of the logarithmic function is a 4-dimensional dot product kap" of two 
null vectors and p " . Vector k" points out the direction of propagation of light ray, while the null vector p" = — z"{s) points out the 
direction of the null characteristic of the gravity field equations. The Lorentz transformation. A" p, from one frame to another changes the 
null vector k" = A" pkl^ , but in order to preserve the Lorentz-invariance of the gravitational time delay At, the null vector p" directed along 
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the body's gravity field must change accordingly p" = A" ^p^ , so that the dot product kap"' = k^^p" remains the same. Hence, not only 
the light undergoes aberration, when one goes from one frame to another, but the null characteristics of the gravitational field in the time 
delay At must change too in the same proportion, if general relativity is valid. In other words, equation ilOl is not the ordinary equation of 
the aberration of light in flat space-time (without gravity field) but a more profound relationship for a curved space-time showing that even in 
the presence of the gravitational field of the moving body, affecting the light propagation, the aberration of light equation remains the same 
as in the flat space-time. This can be true if and only if both the gravitational field perturbation /j^jg and the affine connection F^^ remain 
invariant under the Lorentz group transformation, which is parameterized with the same fundamental speed c as the Lorentz group of the 
underlying electromagnetic wave used in the ranging time-delay experiment. This interpretation is further discussed in more detail elsewhere 
I Kopeikin & Fomalonlll 200611 2007l : [KoDeikin & MakarovlboOTh . 



7 COUPLING OF THE PPN PARAMETERS WITH THE VELOCITY-DEPENDENT TERMS 



7.1 Explicit Coupling 

Equation ([64} describes the Lorentz transformation of the (static) Shapiro time delay from the rest frame of the massive body (Sun, planet) 
to the frame of reference in which the data processing is performed. For we have restricted ourselves with the post-Newtonian expansion of 
the linearized time delay up to the terms which are quadratic with respect to velocity of the moving gravitating body, equation ( I64t can be 
superimposed with the static terms of the second order with respect to the universal gravitational constant G entering equation This is 
because these terms have the same order of magnitude so that we do not need to develop the Lorentz invariant expression for the terms which 
are quadratic with respect to G. We shall also neglect for simplicity the terms which are products of /3- with the PPN parameter fppN because 
Yppn has been already limited by the solar system experiments up to the value not exceeding 10^**. Thus, the product fppNp^ exceeds the 
accuracy of the post-post-Newtonian approximation. 

Our calculation yields the following, Lorentz-invariant equation for the post-post-Newtonian time delay 



At 



+ 



+ 



1 + 



1 + 



Yppn 



k/3- 



Yppn 



kf3- 



2GM 



In 



r\ +r2 + ri2 



2 11) \r\^r2-r\2, 

Yppn \ GM ri2 (ni X /3)-ri + (na x /3)^r2 - (fc x (if{n + ^2) 



(71) 



G^M^ ri2 
£■5 rir2 



c3 rir2 
15 



2Yppn - Pp 



1 + ni • n2 

arccos(ni • n2) 
\ni xn2| 



{2 + Yppn) 



1 +ni •712 



O 



GM 



/3- 



One can immediately observe that the PPN parameter Yppn couples with the velocity terms in front of the logarithmic term. This means that 
the amplitude of the Shapiro delay is effectively sensitive to the linear combination 

f = Yppn - 2fe - IYppnPr + Pr+Pt, (72) 

that will be measured in high-precision space-based experiments like BepiColombo, ASTROD, LATOR, etc. Here and elsewhere, we denote 
respectively jSr = fe • /3 - the radial velocity, and jip = \kx /3\- the transverse velocity of the massive body that deflects the light ray. 

Equation \71i elucidates that the measured value F of the parameter fpPN is affected by the velocity terms, which explicitly present in 
the post-Newtonian expansion of the Shapiro time delay. In case of the ranging gravitational experiment in the field of Sun with the light 
ray grazing the solar limb, one has d = Rfi, = 7 x lO'" cm - the solar radius, and ;'g = 3 x 10^ cm - the Schwarzschild radius of the Sun. 
The Sun, in moving in its orbit around the barycenter, has an average distance of 1.1 Rq from it but may be as far as 2.3 Rq. The orbital 
path of the Sun about the barycenter traces out a curve that is closely resemble an epitrochoid - three-lobed rosette, with three large and 
three small loops - with a loop period of 9 to 14 y ears. Fifteen successive orbits comprise a 179-year cycle of the solar motion around the 
barycenter 1 Fairbridge & Shirlevlll987 ; Jose 1965 ) - the duration, which is also the time taken for the planets to occupy approximately the 
same positions again relative to each other and the Sun. The solar velocity vq with respect to the barycenter of the solar system can reach 
maximal value of 15.8 m/s giving rise to jS© = Vq Ic = 5.3 x 10~^. Because space mi ssio ns LA TOR and ASTROD are going to measure fpPN 
parameter with a precision approaching to 10"' dTurvshev. Shao & Nordtved^l 20041 : iNilboOTh . the explicit velocity-dependent correction to 
the Shapiro time delay in the solar gravitational field must be apparently taken into account. Current indeterminacy in the solar velocity 
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vector is about 0.366 m/day Jpitievalboosh that yields an error of Afe ~ 1.4 X 10^14 This error is comparable with the contribution of the 



second-order velocity terms jSj, ^ 2.8 x 10^'^. However, they are too small and can be neglected in the measurement of yppf^. 

Coupling of the velocity-dependent terms with parameters ^ppN and dppN can be understood after making expansion of high-order 
terms in equation l l71b with respect to the impact parameter of the light ray rf = |fcxri| = |fexri| that is assumed to be small: d <tiri, 
d <^r2. The unit vectors ni and n2 can be decomposed in the post-post-Newtonian terms as follows 



ni = — fccos 6] +nsin 01 , 
712 = k cos 62 + n sin 62 , 

where the unit vector n is directed from the massive body to the light-ray trajectory along the impact parameter: d - 
introduce the deflection angle Q defined as 

n\ n2 = cos(;r — 0) = — cos Q . 



(73) 
(74) 

: dn. It is convenient to 



(75) 



One can easily observe that 6 = 0] + 02- Practically all gravitational ranging experiments are done in the small-angle approximation, when 
<C 1, 9i <C 1, 02 ^ 1- In this approximation, one has 



q1 , 
l+ni-n2 = — + 0\e' 

(ni x/3)Vi+(n2X/3)2r2-(fcx/3)2(r,+r2) = 0d(^j5~-j5^ 
Substituting equations i72\ . ll75t-(l77b to equation (inj yields 

ri+r2 + ri2 



6 



At 



ri + ?-2 - ri2 



(76) 
(77) 

(78) 



+ 



G^M^ ri2 
ri 1-2 



15 



2fppN-j3, 



>PPN- 



- SppN 



2(2 + Yppn)- 



02 



O 



GM , 



where we have introduced a new notation 
fppN\ 16 d 



SppN = Sj 



'PPN - 



J 3k] 



(79) 



and denoted Vg = 2GM/c~ - the Schwarzschild radius of the massive body deflecting the light ray. Explicit contribution of the solar velocity 
terms to the parameter 8ppM ca n achieve 1.1 x 10~^ that is much less than the precision of measurement of the PPN parameter Spp^j in 
LATOR and ASTROD missions 1 Plowman & HeflingJlOOeh and can be currently neglected. 

We recall to the reader that in scalar-tensor theory o f gravity parameter ^ppN can not be determined separately from Sppf^ as they appear 
in the linear combination —pppN + 3/4dppM- Following 1 Plowman & Hellingsl 20061) we assume that Pppn is determined from other kind of 
gravitational experiments, and eliminate it from the fitting procedure. 



7.2 Implicit Coupling 

In the previous section we have made an explicit post-Newtonian expansion of the ranging time delay in powers of the velocity-tracking 
parameter j6 = v/c. This post-Newtonian expansion is shown in equation J7U . It looks like the only place, where the linear velocity correction 
to the Shapiro delay appears, is in front of the logarithmic term. However, a scrutiny analysis reveals that the linear velocity-dependent 
correction is also present implicitly in the argument of the logarithmic function. Indeed, distances ri = |xi — z{ti)\ and r2 = \xi — z{t2)\ 
depend on two positions of the massive body taken at two different instants of time, ti and t2- The body moves as light propagates from the 
point of emission xi to the point of observation X2, so that the coordinates of the body are not arbitrary but connected through a relationship 

z{t2) = z{ti)+v{t2-ti) , (80) 

which, indeed, shows that the velocity of the body is involved in calculation of the numerical value of the argument of the time-delay 
logarithm. 
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Though this dependence on the velocity of the massive body is implicit, it definitely affects the measured values of the PPN parameters 
and makes their values biased in case if either general relativity is invali d or if the numerical code used for data processing of the ranging 
experiment, does not incorporate the solar system ephemeris properly jKopeikin et al.ll 20071) . Let us show how this impact on the PPN 
parameters can happen. 

To this end we shall assume that the light ray passes at a minimal distance d from the body at the time of the closest approach t^, which 
is defined in the ap proximation of the unpertu rbed light-ray trajectory, x{t) = xi + k(t — 1\) for t ^ t\ oi x(t) = X2 + k{t — 12) for t ^ t2, 
from the condition 1 Klioner & Kopeikin|[l992h 



dt 



-0, 



(81) 



where x{t) = xi+k(t — ti) is the (unperturbed) light-ray trajectory, and z(t) = z{ti) + v{t — ti) is the body's world line in the approximation 
of a straight line, uniform motion. Taking the time derivative and solving the equation yield 



(T ■ ri 



--t2- 



(T ■ r2 



(82) 



c\k-f3\ c\k-f3\ ' 

where the unit vector ct has been defined in equation i57i . The post-Newtonian expansion of various distances near the time of the closest 
approach gives us 



r2* 



l-{/3-nu) — + — 

1- /3-n2*) h — 

ri* 2 



l-/3fc 



(/3 X kf 



(83) 
(84) 
(85) 



where / j = c(ti—t^), I2 = c{t2 — t^), the unit vectors ni^ = ri^/ri^, n2, = r2»/r2,, and distances ri^ = xi — z{t^), r2» = X2 — z(t^). 

We substitute now the post-Newtonian expansions ll83t-J85t to the logarithmic function of the Shapiro time delay and apply the small- 
angle approximation. It will yield 



In 



n+r2 + ri2 
n +r2~r\2 



■ in 



rit + r2t+r\ Ird^ k (3 



[l + O(j8) + O(0,) 



(86) 



where 6* is the angle between two vectors ni^ and n2» defined as ni* ■ n2» = cos(;r — 6,). 

The post-Newtonian expansion of the ranging delay in the vicinity of the time of the closest approach of the light ray to the massive body 
reveals that the parameter SppN is affected by the first-order velocity terms from equation ( I86t . Specifically, taking into account equation 
l l86t allows us to write down the ranging delay in the following form 



At 



+ 



where 



(2+r 



GM 



In 



ri*r2* 



ru +r2* —r 
15 
~4 



(87) 



'27PPN - fipPN + ^ 



n 

0^ 



2{2 + fppN 



o 



(GM 



A =8, 



'PPN- 



2 I in Tg 



(88) 



The last term in equation l l88t can amount to 0.02, which exceeds the expected accu racy of measuring the PPN parameter 8ppf^ with 
LATOR/ASTROD missions by a factor of 10 as follows from l lPlowman & Hellingsl200i) . This clearly indicates the necessity of inclusion of 
the velocity-dependent post-Newtonian corrections to the data analysis of the high-precise time delay and ranging gravitational experiments. 



© 2009 RAS, MNRAS OOO.fTlfTsl 



14 S. M. Kopeikin 



8 RANGING EXPERIMENTS AND LORENTZ INVARIANCE OF GRAVITY. 



In special relativity, where the Minkowski geometry represents a flat space-time, the Lorentz symmetry is a global symmetry consisting of 
rotations and boosts. However, in curved space-time, in the most general case, the Lorentz symmetry is a local symmetry that transforms local 
vectors and tensors in the tangent (co-tangent) space at each space-time point. Nonetheless, general rela tivity admits the global Lorentz sym- 
metry, at least, for isolated astronomical systems residing in asymptotically-flat space-time i Focklll964 ). This asymptotic Lorentz symmetry 
of gravitational field can be traced in the invariant nature of the gravitational Lienard-Wiechert potentials given by equation JlOb . which are 
solutions of the linearized Einstein equations. The asymptotic Minkowskian space-time for isolated systems defines the background manifold 
for gravitational field perturbations, h^j^a , and must have the same null-cone structure as the local tange nt space-time, which is defined by mo- 



tion o f light particles (photons). However, this theoretical argument is a matter of experimental study jKopeikiijboOll : iFomalont & Kopeikinl 
2003h . 

Ranging time-delay experiments are, perhaps, the best experimental technique for making such test. This is because they operate with 
the gauge-invariant fundamental field of the Maxwell theory having well-established and unambiguous physical properties. Propagation of 
radio (light) signals traces the local structure of the null cone hypersurface all the way from the point of emission down to the point of its 
observation. Now, if the massive body, which deflects radio (light) signals, is static with respect to observer, one can not draw any conclusion 
on the asymptotic structure of the space-time manifold and on whether its Lorentz symmetry is compatible with the Lorentz symmetry of the 
light cone. This is because the gravitational interaction of the body with the radio (light) signal is realized in the form of the instantaneous 
Coulomb-like gravitational force with having no time derivatives of the gravitational potentials been involved. However, if the massive body 
is moving with respect to observer as light propagates, its gravitational force is not instantaneous and must propagate on the hypersurface of 
the null cone of the asymptotic Minkowskian space-time as it is described by the Lienard-Wiechert gravitational potentials l IlOl l. The terms 
in the ranging time-delay l l32t depending on both the translational velocity (3 = v /c of the massive body and the retarded time s, originate 
from the time derivatives of the gravitational potentials and characterize the global Lorentz sym metry of the gravitational field . Therefore, 
measurement of these terms in the ranging time-delay experiments has a fundamental significance i Kopeikin & Fomalonlll2006l) . 

Currently, there is a growing int erest of theore tical physicists to gravitational theories where the global Lorentz symmetry of gravitational 
field can be spontaneously violated jBluhmll200^ . This is motivated by the need of unification of the gravity field with other fundamental 
interactions. These theories introduce additio nal long-range fields to the gravitational Lagrangian, which des troy the symmetry between 
the, so-called, observer and particle invariance I Kosteleckv & Pottingll995 : Colladav & Kosteleckvlll99~ . 1998 ). Interaction terms involving 
these fields appear also in the equations of motion of test particles. It is the interaction with these fields that can lead to physical effects of 
the broken Lorentz symmetry that can be tested in experiments. Outcome of these experiments depends crucially on the assumptions made 
about the structure of the additional terms in the gravitational Lagrangian and the numerical value of the coupling constants of these fields 
with matter. On the other hand, the measurement of the post-Newtonian velocity-dependent and/or retarded-time corrections in the ranging 
time-delay experiments does not depend on any additional assumptions and relies solely on general relativistic prediction of how the radio 
(light) signals propagate in time-dependent gravitational fields. 

It is remarkable that current technology already allows us to measure the velocity-dependent and/or retarded-time post-Newtonian 
corrections in the ra nging time-delay experiments conducted in the solar system. The most notable experiment had been done in 2002 with 
the VLBI technique jFomalont & Kopeikinll2003r) . It measured the retard ed component o f the n e ar-zone gravitational fie ld of Jupiter via its 
impact on the magnitude of the deflection angle of light from a quasar i Kopeikinl I2OO1L l2004h . iFomalont et all hoogh have repeated this 
retardation of gravity experi ment in 2009 by making use of the close encounter of Jupiter and S aturn with quasars in the plane of the sky. 

The Cassini experiment i Bertotti. less & TortoralfioOS : Anderson. Lau & GiampierilEooi) is also sensitive to t he time-dependent per- 
turbation of gravitational field of the Sun c aused by its orbital motion around the barycenter of the solar system jKopeikin et al.ll 2007 : 
Bertotti. Ashbv & Iessl20oijKopeikinl2009l) . However, its detection requires re-processing of the Cas sini data in order to separate the Cassin i 
measurement of PPN parameter y/jpA? from the gravitomagnetic deflection of light by the moving Sun jKopeikin et alj 2007l : lKopeikinll2009l) . 



9 RANGING DELAY IN THE NASA ORBIT DETERMINATION PROGRAM 



Relativistic ranging time delay, incorporated to the NASA ODP code, was originally calculated by iMoven ( 120030 under assumption that 
the gravitating body that deflects light, does not move. Regarding the Sun, it means that the ODP code derives the ranging delay in the 
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heliocentric frame. Let us introduce the heliocentric coordinates X" = {cT,X'), and use notation x" = (ct,x') for the barycentric coordinates 
of the solar system, which origin is at the center of mass of the solar system. The Sun moves with respect to the barycentric frame with 
velocity vq = dxQ /dt amounting t o ~ 15 m/s. Though th is velocity looks small, it can not be neglected in such high-precision relativity 
experiments as, for example, Cassini i Kopeikin et al.l 2007 ). A legitimate question arises whether the ODP code accounts for the solar motion 
or not. We demonstrate in this appendix that the ranging time delay in the ODP code is consistent with general relativity in the linear-velocity 
approximation, but it fails to take into account the quadratic velocity terms properly. Thus, more advanced theoretical development of the 
ODP code is required. 

The ranging time delay in the heliocentric coordinates with the Sun located at the origin of this frame, follows directly from equation 
( 1441 ) after making use of the heliocentric coordinates. It reads 



AT 



^\X2-Xi\+AT 

2GMf:, , 

■In 



,.3 



R2+R1 -Ri 



(89) 
(90) 



where X2 and are the heliocentric coordinates of observer and emitter respectively, distance of the emitter from the Sun is R2 = |-X^2|> 
distance of the observer from the Sun is /?i = Xi, and 12 = | X2 — Xi \ is the null heliocentric distance between the emitter and observer. 
This equation coincides exactly (after reconc iling our and Moyer's notati ons for distances) with the ODP time-delay equation (8-38) given 
in section 8 of the ODP manual 1 Move J 20031) on page 8-19. Movei 1 2003 ) had transformed the argument of the logarithm in the heliocentric 
ranging delay l |90l > to the barycentric frame by making use of substitutions 



X2 r2 = a;2 — a;0(f2) , Xi ^ r2 ~ xi — X(^,{t\) . (91) 

The ODP manual ( lMovej|2003l) does not provide any evidence that these substitutions in the ranging time delay ( |90l > are consistent with gen- 
eral relativity and do not violate the Lorentz symmetry. Nonetheless, comparison of equations ( 1901 ), l |91t with the post-Newtonian expression 
i64i for the ranging delay demonstrates that equations ( 1911 ) are legitimate transformations from the heliocentric to the barycentric frame in 
the sense that they take into account velocity of the Sun in the ranging time delay in the linearized, post-Newtonian term following the static 
Shapiro time delay. 

Equation ( 1641 ) also shows that the ODP code is missing the velocity-dependent term in front of the logarithmic function in equation ( |90t . 
The ranging time delay in the heliocentric and barycentric frames must be related by the simple equation 

At={l-k-l3e)AT , (92) 

which is a linearized version of equation i64\ that was derived by Kopeikin & Schafei ( 19991) . We conclude that the ODP code used by NASA 
for navigation of spacecrafts in deep space, is missing a high-order velocity-dependent corrections to the Shapiro time delay and can not be 
used for processing and unambiguous interpretation of near-future ranging experiments in the solar system. A corresponding relativistic 
modification and re-parametrization of the ODP code based on equations of the present paper is highly required. 

Equation J92I ) has been also derived by Bertotti. Ashbv & less 1 200 sh who claimed that the velocity-dependent terms appear in the time 
delay only in front of the logarithmic function in equation ( 1921 ). As we have shown in section [7^ the argument of the logarithm in equation 
( 1631 ) also contains terms depending on velocity v of the gravitating body, which are implicitly present in the definition of the distance ri2. 
This distance is calculated between two spatial points separated by the time interval required by light to travel between the point of emission 
and observation respectively (see equations (56)-(58)). Coordinates z{ti) and z(t2) are not the same because the gravitating body is moving. 
These coordinates are related by means of the equation (72), which demonstrates that velocity v of the gravitating body must be known in 
order to calculate the distance r^. Because one has to rely upon equation (72) in the ODP data processing algorithm, the post-Newtonian 
expansion of distance ri2 yields 

n2 = r-r-f3 + 0{l3-) , (93) 

where the null distance r = \r\ is defined in equation i35\ . It follows that the distances ri2 and r entering equation i93\ are not the same 
quantities as they differ by terms of the order of v/c. Equation \93\ reduces the ranging delay i92\ to the following form 

r — rj3 



At 



(l-fc./3)^ln 



fl + ri 



r2 + n -r + r-/3 



(94) 
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which has been derived in our paper i Kopeikin et al.ll 20 07). Bertotti. Ashbv & les^ 1 20081) claimed that the expression l lgit for the ranging 



time d elay does not appear in the ODP manual (|Moyejj|20oJ and is not allowed for theoretical analysis of the Cassini experiment as we 
did in jKopeikin et al.ll 20071) . However, expression l |94| l is exactly the same function AT" given in the ODP manual but expressed, instead of 
distance r\2, in terms of the distance r and velocity of the Sun, v, via self-consistent mathematical transformation J93b . For this reason, the 
two express ions are mathem atically equivalent and either of them can be used in data processing of the ranging observations of the Cassini 
experiment jKopeikiijEooi) . 
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Light-ray deflecting body 



Figure 1. Ranging time delay experiment. Electromagnetic signal is emitted at distance r\ from the massive body, passes by it at the minimal distance d, and 
is received by observer at distance ri. The emitter, observer, and the massive body move with respect to each other as the electromagnetic signal propagates. 
This makes the ranging delay experiment sensitive to the null cone structure of space-time in general relativity and modifies the Shapiro time delay. 
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